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Dephasing and Hyperfine Interaction in Carbon Nanotubes 
Double Quantum Dots: The Clean Limit 
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We consider theoretically ^^C-hyperfine interaction induced dephasing in carbon nanotubes double 
quantum dots with curvature induced spin-orbit coupling. For two electrons initially occupying a 
single dot, we calculate the average return probability after separation into the two dots, which 
have random nuclear-spin configurations. We focus on the long time saturation value of the return 
probability. Poo. Because of the valley degree of freedom, the analysis is more complex than in, for 
example, GaAs quantum dots, which have two distinct Poo values depending on the magnetic field. 
Here the prepared state and the measured state is non-unique because two electrons in the same dot 
are allowed in six different states. Moreover, for one electron in each dot sixteen states exist and 
therefore are available for being mixed by the hyperfine field. The return probability experiment is 
found to be strongly dependent on the prepared state, on the external magnetic field — both Zeeman 
and orbital effects — and on the spin-orbit splitting. The lowest saturation value, being Poo = 1/3, 
occurs at zero magnetic field for nanotubes with spin-orbit coupling and the initial state being the 
groundstate, this situation is equivalent to double dots without the valley degree of freedom. In 
total, we report nine dynamically different situations that give Poo = 1/3, 3/8, 2/5, 1/2 and for 
valley anti-symmetric prepared states in an axial magnetic field, Poo = 1. When the groundstate is 
prepared the ratio between the spin-orbit splitting and the Zeeman energy due to a perpendicular 
magnetic field can tune the effective hyperfine field continuously from being three dimensional to 
two dimensional giving saturation values from Poo = 1/3 to 3/8. 

PACS numbers: 85.35.Kt, 81.05. ue, 73.21.La, 31.30.Gs, 03.65.Yz 



I. INTRODUCTION 



Quantum dots are attractive candidates for imple- 
menting qubits; in particular carbon-based materials 
such as nanotubes (CNTs) and graphene provide the ad- 
vantage of a weaker hyperfine interaction. In general the 
samples have high concentration of ^^C isotopes (spin 
zero) and low concentration of ^"^C isotopes (spin-1/2), 
and thus the coupling between the confined electron spin 
and the nuclei spins is small. However, the valley de- 
gree of freedom increases the number of available few- 
particle states in the quantum dots, leading to a more 
complex system, as compared to, e.g., GaAs quantum 
dots. Breaking of this 4-fold spin and valley degeneracy 
has been predicted for nanotubes due to a curvature- 
induced spin-orbit coupling,^"'' which has been confirmed 
experimentally.^^^ The spin-orbit split spectrum leads to 
doubly degenerate states that can be used as qubit states. 
For example, by taking advantage of the strong diamag- 
netic effects in an axial magnetic field, a spin qubit at 
larger magnetic fields can be defined. * Alternatively, the 
low field spin-orbit entangled Kramers pair can be used 
as the qubit, which has been proposed as an electrically 
manipulatable qubit in bend nanotubes.^ 

As pointed out by Loss and DiVincenzo,^^ one of 
the requirements for quantum dot-based quantum com- 
putation is the ability to control the exchange inter- 
action, which has been successfully demonstrated in 
GaAs double quantum dot (DQD) setups. ^^'^^ Simi- 



lar devices have more recently been fabricated in CNT 



systems 



6,13 



which have also attracted considerable the- 



oretically attention. 



14-17 



Another key requirement is long decoherence times. 
Experimental determination of the spin decoherence 
times in DQD is done by converting the spin information 
into charge information, which can be measured using 
on-chip charge detectors. The charge state of the DQD 
is characterized by the number of electrons in the left and 
in the right dot, (N^jNn), which are controlled by gate 
voltages applied to the left and to the right dot, Vl and 
Vff, producing a map of the equilibrium charge states of 
the double quantum dot, as illustrated in Fig. 1(a). The 
figure also shows the cycle that allows measurement of 
the dephasing time. This is done by initially preparing 
the system in a state with two electrons in one dot, then 
separating them for a certain time with one electron in 
each dot — while the initial state can be affected by the 
environment — and finally by measuring the probability 
for the electrons to return to the original dot. 

In this paper, we study the return probability experi- 
ment, which measures the characteristic decay time, T^ , 
for the return probability averaged over many cycles, as 
well as the saturation value, Poo, see Fig. 1(b). In the case 
when the dominant time evolution is due to hyperfine in- 
teraction with a nuclear spin background that changes 
between cycles, the r| time is a measure of the inho- 
mogeneous broadening of the spin state. This is well- 
studied in two-dimensional electron gases defined quan- 
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FIG. 1. Return probability experiment, (a) Cycle applied in 
gate voltage space {Vl,Vr) in each single-shot measurement, 
(b) Sketch of the return probability, P{rs), as a function of 
the time waited, Ts, at "e"; the saturation value. Poo, and 
the decaying time, r^, are included, (c) Spectrum of the (1,1) 
and (0,2) two-particle states (taken for Ago ^ 0, tunneling 
gap. At, hyperfine characteristic energy, ah, and B = 0) as a 
function of the detuning, e; arrows and detuning values, £p, £o 
and £m refer to the single-shot cycle's stages presented in (a). 
We label with squares ( "e.g" ) the states that the system may 
occupy at the stages, "p","e" and "m" if the (0,2) highest 
energy excited state is prepared; the possible outcomes at Sm 
would be different if the prepared state, |^p), had been one 
of the other five available (0,2) states. 

turn dots."'^^"^^ 

Here, due to the valley degree of freedom, the prepared 
state with two electrons in one dot is not restricted to 
one possibility, as the spin singlet in a GaAs dot, but 
six states. Furthermore, once the electrons are sepa- 
rated, there are sixteen available states instead of four 
as in a GaAs double dot (the spin singlet and the spin 
triplets). As mentioned above, the spin-orbit coupling 
and the diamagnetic effect of an external magnetic field 
split up these manifolds of states. In Fig. 1(c) we present, 
as a function of the detuning, e, the spectrum of the 
six (0,2) states, the sixteen (1,1) states and the mixing 
between them due to single-particle tunneling. The re- 
lation between the detuning and the gate voltages, the 
nature of the mixing of the states, and the description 
of the stages in the single-shot measurement cycle, are 
presented in detail in Sec. II and Sec. III. 

So far one experimental study of Tj* nieasurements has 
been reported^ on samples with high concentration of 
^^C. The result shows a saturation value in the return 
probability of ~ 0.17, which is not understood. The the- 
ory for non-valley degenerated DQDs predicts (as lower 
bounds) 1/3 and 1/2 in the low and high magnetic field 
regimes, respectively. A multivalley case was recently in- 
vestigated for a silicon double dots.^^ Depending on the 
prepared state, either a GaAs-like behavior or hyperfine 
immune states may be found. The hyperfine interaction 



considered did not, however, involve valley mixing terms 
in contrast to the C-based dots we treat here. 

The motivation of our detailed study is to predict 
the expected return probabilities for carbon based sys- 
tems with hyperfine coupling to spin 1/2 ^■^C nuclei. In 
these graphene-based systems, in contrast to Si double 
dots, one must take into account that the hyperfine in- 
teraction affects both the spin and the valley degrees of 
freedom. ^^'^^ We include those hyperfine valley mixing 
effects but we do not include disorder induced spin con- 
serving valley mixing, which is presented in a separate 
publication.^^ Furthermore, we include direct Coulomb 
interaction, but not Coulomb exchange which is expected 
to be a small effect. ^^ 

We also work in the limit of large detuning so that the 
tunneling exchange on the (1,1) states is much smaller 
than the hyperfine field characteristic energy and there- 
fore we obtain the lower bounds for Poo- On the other 
hand, if the tunneling exchange is important, the de- 
generacies for zero hyperfine are reduced diminishing the 
effectiveness of the hyperfine-induced mixing and there- 
fore also increasing the saturation values of the return 
probability. In those situations Poo grows continuously, 
as a function of the tunneling exchange, from the zero- 
exchange value up to one.^" 

We start from a simple model for an isolated quantum 
dot and construct from this the two-electron wave func- 
tions. We study the cases of large and small spin-orbit 
coupling. The result is found to be dependent on the pre- 
pared state and on the external magnetic field. Notably, 
for some situations the two-electron wave function is al- 
most not dephased by the hyperfine field. We show that 
besides the usual saturation values of the return probabil- 
ity, 1/3 and 1/2 — well known in DQDs without the valley 
degree of freedom and for zero tunneling exchange — other 
values can be observed; namely, 3/8, 0.4 and 1. In addi- 
tion, for nanotubes with spin-orbit coupling, an applied 
magnetic field in a direction perpendicular to the tube 
axis can tune the saturation value between 1/3 and 3/8. 

The paper is organized as follows. In section II, we 
describe the double dot model and the four special cases 
used in the paper. Section III describes the experiment 
and the methods used to calculate the return probabil- 
ities, with results presented in Section IV. Finally, con- 
clusions and summary are found in Section V. 



II. QUANTUM DOT MODEL 

A. Single dot 

We consider semiconducting tubes where the bandgap 
is due to either chirality or, for nominally metallic tubes, 
to curvature. ^^'^^ The semiconducting properties allow 
electrons to be confined in a gate-defined potential.^' ^''"^^ 
This potential is assumed smooth on the scale of the in- 
teratomic distance, conserving the valley index. Both 
single* and double quantum dots^"*"^^ have been studied 



in this approximation. An important effect of curvature 
in nanotubes is that it leads to spin-orbit interaction that 
couples the valley index with the spin in the longitudinal 
direction. "'^"^ 

The Hamiltonian describing the spin and valley degrees 
of freedom, excluding the hyperfine interaction, reads 



-ffsv — ^S( 



H, + H, 



orb; 



(1) 



where Hgo is the spin-orbit coupling term, Hg is the Zee- 
man interaction and iJorb is the diamagnetic effect of the 
magnetic field. For the specific cases considered below, 
not all terms in Eq.(l) are present. 
The spin-orbit coupling term is 



ffs, 



1 



DT3CT11 



(2) 



where Ago is the spin-orbit energy splitting, a\\ is the 
spin operator along the direction of the tube axis and 
(7a;, ay, az ((To) and ti, T2, T3 (tq) are the Pauli (identity) 
matrices in spin and valley space, respectively. For the 
valley degree of freedom, we use r = K, K' to identify 
the -1-1 and —1 eigenstates of T3. In the following, unless 
otherwise stated, the spin quantization axis — a =t, i or 
equivalently, a = ± — is taken along the direction of the 
total magnetic field. The value of the spin-orbit split- 
ting depends on the nanotube's chiral vector and on the 
electron filling.^ 

The Zeeman energy due to an external magnetic field 



B, 



[BxiBy^Bz) is 



Hs = TidsfJ-BTo (Bex ■ cr) , 



(3) 



where gs=2 is the usual gyromagnetic factor. The com- 
ponent of the magnetic field parallel to the tube axis, B\\ , 
gives rise to a strong diamagnetic effect: 



Hn 



goiht^BBnT^ao, 



(4) 



where the orbital g-factor, goih, depends on the size of the 
nanotube'^^ and it is bigger than gs/2 in the typical case: 
nanotubes with radius greater than 1 nanometer. ^^'' To 
simplify the notation, we define the following two energy 
scales: 



Es — .gs Mb I Bex I, 

Eoih — 25orbM-B|-B|||. 

1. Four special cases 



(5) 
(6) 



Here we identify four special cases of single quantum 
dots, representing different physical realizations, see also 
Fig. 2. The important classification is the splitting of 
the dot state compared with the hyperfine energy. The 
relevant energy scales are Zeeman splitting, Eg, orbital 
splitting, i^orb, and spin-orbit splitting. Ago; we only con- 
sider each one of them when they become much bigger 
than the hyperfine interaction. Based on this, we define: 
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FIG. 2. Single-particle energies as a function of magnetic field 
for the quantum dots of cases A, B, C and D presented in the 
text. For A, B and C the eigenstates are |ro-) because the spin 
projection along the direction of the external magnetic field 
is a good quantum number, the corresponding eigenenergies 
are given in Eq.(lO). The parameters are: (a) Aao = i?orb=0 
and Ss ^ 0; (b) A,o = 0, Sorb ^ and E^ ^ 0; (c) A,o / 0, 
Sorb / 0, Ss / and B_l = 0. (d) As in (c), A^o 7^ 0, but the 
external magnetic field is perpendicular to the tube's axis, so 
-Bs 7^ and JSorb = 0; the eigenfunctions, |t, ±), change their 
spin projection according to the ratio between the Zeeman 
energy and the spin-orbit splitting, see details in Sec. IV D. 



Case A) No spin-orbit coupling and no orbital mag- 
netism, only the Zeeman energy E^ is considered. This 
is relevant for nanotube dots with small spin-orbit en- 
ergy in a perpendicular magnetic field (so that ii^orb — 0) 
or, alternatively, graphene dots in an in-plane magnetic 
field. "^^ The single-particle spectrum as a function of the 
magnetic field is shown in Fig. 2(a). Results for valley 
degenerated quantum dots are presented in Sec. IV A. 

Case B) Quantum dot with no spin-orbit coupling 
and finite orbital magnetism; this is the situation for 
nanotubes with small spin-orbit splitting in a magnetic 
field with a parallel component. In this case, the ener- 
gies Es and -Eorb are both relevant. Figure 2(b) shows 



the spectrum as a function of the magnetic field, which 
breaks the 4-fold degeneracy. Since we have assumed that 
-Sorb > Eg, which is likely the case if the total magnetic 
field is parallel to the tube axis,^"^ the two highest en- 
ergy solutions have r = K, whereas the two lowest energy 
solutions have t = K'. Results for the current case are 
presented in Sec. IV B. 

Case C) Quantum dot in a nanotube with spin-orbit 
coupling and external magnetic field parallel to the tube 
axis; all three energy scales Ago, E^ and ii^orb are rele- 
vant. At zero magnetic field, in contrast to cases A and 
B, the spin-orbit coupling breaks the four- fold degener- 
acy, which results in two Kramers doublet. The energy of 
the two Kramers doublets are ±Aso/2, as shown in 2(c), 
which depicts the spectrum as a function of the magnetic 
field. For two finite fields some degeneracies are recov- 
ered. Results for this situation are presented in Sec. IV C. 

Case D) Here the system has a finite spin-orbit cou- 
pling and the magnetic field is perpendicular to the tube 
axis, therefore, diamagnetic effects are absent. Even for 
finite Zeeman energy. Eg, the doublets are not split be- 
cause the magnetic field cannot couple opposite valley, 
see Fig. 2(d). The ratio between E^ and Ago controls the 
spin projection of the solutions; the bigger the magnetic 
field the more similar to case A the solutions become be- 
cause the spin of the eigenstates tend to align with the 
magnetic field. This quantum dot solutions and the re- 
turn probability results are presented in Sec. IV D. 



The double dot 



The double dot single-particle Hamiltonian is 



-"dqd 



<R 



H'^, 



where ^l^^j^ includes the Hamiltonians for the two iso- 
lated dots and Hrj^ is the single-particle tunneling be- 
tween the two dots. In what follows we use the su- 
perindex Ip or 2p to distinguish single-particle and two- 
particle operators. We introduce the Pauli (identity) ma- 
trices i,i (^o) in left-right space, with ±1 eigenvalues of 
^3 for L/R, respectively. In this notation, the tunneling 
part of the Hamiltonian becomes 



H^ = -tS,lTQ(TQ. 



(8) 



This inter-dot tunneling is assumed to be valley and spin 
conserving, because the gate voltage defined confining 
potential is non-magnetic and it is assumed smooth on 
the lattice scale. We also assume that the tunneling am- 
plitude does not depend on the quantum numbers r and 
(T, which is valid as long as the height of the potential 
barrier is much bigger than the detuning, the spin-orbit 
and magnetic splittings. ^^■"'^^ 

The isolated left dot plus right dot single-particle 
Hamiltonian is 



rrlp 

-"l,r 



ToO-Q- 



Co(eL + eR) +6(eL - er) 



+ Co^s' 



(9) 



where the last term is the valley and spin Hamiltonian in 
Eq. (1) which is identical in the two quantum dots. The 
effects of the gate voltages are introduced as the energy 
shifts eL and er for the left and the right dot, respectively. 



C. Two-particle basis states, no tunneling 

A single-particle basis set can be generated by the eight 
states I^Tf), with ^ = L,R, r — K, K' and the spin 
projection a —^, \, is taken along the direction of the 
magnetic field. If Ago — 0, or if the magnetic field is 
parallel to the tube's axis, \iTa) are eigenstates of the 
single-particle Hamiltonian H^^^ with eigenenergies 



E^ra = {EsCr + EorhT - AgoTCr) /2. 



(10) 



Using these states we build two-particle Slater deter- 
minants with quantum numbers ^ra and ^'t' a' as follows 



?'t'<t'\ _ 



1 

71 



(|eTa),|rrV),-|C'rV),|Cra),) 



(11) 

Since the single-particle basis has eight elements, the two- 
particle basis has 28 states (28 = 8!/(6!2!)). However, for 
the return probability due to energetic reasons we do not 
include the states with two electrons in the left dot, i.e., 
the (2,0) charge configuration, which leaves 22 states. 

In general, the single-particle eigenstates can differ 
from I^Tcr); we label the single-particle eigenstates of iJgv 
as |n), with energies En using the index n = 1,...,4; 
then, the six (0,2) eigenstates and their eigenenergies are 



(7) (0,2) 



Rn' 
Rn 



E, 



(0,2) _ 



2eR + ^„ + £;„. +C/RR, (12) 



where n, n' = 1,...,4 and n < n' . Here we included 
electron-electron interaction, represented by the right dot 
charging energy Urr. We have not included Coulomb 
exchange since it is expected to be small. "^'^^ With one 
electron in each dot the sixteen (1,1) eigenstates and their 
corresponding eigenenergies are 



(1,1) 



Ln' 
Rn 



E, 



(14) 



e^+e-R+En+En'+Ui^R, (13) 



where J/lr is the inter-dot Coulomb repulsion. 



D. Inter-dot tunneling 

The single-particle Hamiltonian of Eq.(8) preserves the 
spin and valley degrees of freedom in the inter-dot tun- 
neling (L o R) , therefore, the tunneling Hamiltonian can 
be rewritten in terms of the single-particle eigenstates as: 



H]? 



-i^(|Ln)(Rn| + |Rn)(Ln|)^ 



(14) 



When acting with the tunneling Hamiltonian on the (0,2) 
two-particle basis states it gives, 



H; 



2p 



Rn 



= -t 



Rn' 
Ln 



Ln 
Rn 



(15) 



i.e., a combination of (1,1) Slater determinants associated 
with those single-particle states. Thus, any given pair 

n, n! gives a 3 x 3 Hamiltonian. In the basis. 



Rn' 
Rn 



Ln' 
Rn 



and 



L^ ), this Hamiltonian matrix becomes 



Hnn' — EaV + En + En' 



2 

^ 2 

-t 



t -t 





(16) 



where the detuning, e, and the average energy, Eav , have 
been defined as 



Eav 



E. 



(14) 



E. 



(0,2) 



cl 



(1,1) 



E 

1 



K 



(0,2; 



(En 



{el + 3e_R + t/LR + %,r) . 



C/rr + t^LR, (17) 

^En') 

(18) 



For the mixing of the (0,2) and the (1,1) states, the global 
energy shift EAv+En+En' is irrelevant, and therefore we 
choose Eav — from this point on. The Hamiltonian of 
Eq.(16) can in fact be reduced to a 2 by 2 system because 
the following (1,1) combination. 



LR^ 



odd 



1 



Rn' 

Ln 



Ln' 
Rn 



(19) 



p(i,i) 



is an eigenstate with energy independent of t, E^ ^, ^^^ 
^ + En + En' ■ In the Pauli blockade language this is a 
blocked state. The left-right (LR) symmetry of this state 
is evident when writing it as a product in LR- and nn'— 
spaces: 



LR^ 



'odd 



^5 y I ln,n'/ ) 



(20) 



where we introduced the notation of a "singlet" in LR- 
space. It is convenient to introduce also the triplet states 
in this space, and we define 



|SL^) 



|L>,|R)2-|R>i|L). 



I LRv _ |L),|R), + |R),|L), 



(21b) 



^0 

|tlR) = |l),|l),, 
|tl^) = |r>,|r),. 

Similarly, we have defined singlet-likc and triplet-like 
functions in valley and spin spaces (not shown). Finally, 
for the combined quantum number n, we define 
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FIG. 3. Spectrum of the (0,2) and (1,1) states in the presence 
of inter-dot tunneling, t, as a function of the detuning, e. The 
plot focuses on the three states based on the single-particle 
states \n) and \n') for one of the six independent cases (with 
n^n'). The avoided crossing, with gap At = 2f\/2, at £ = 
is a signature of the mixing between (0,2) and (1,1) states. 
The (1,1) odd combination of Slater determinants involves 
a singlet in LR-space (see Eqs.(20) and (24d)) and it is not 
affected by the tunneling. As shown in Eqs.(15) and (26), 
the (0,2) state (i.e., a T^^ state) is mixed only with the even 
combination of the shown (1,1) Slater determinants (i.e., a 
To^ state). 

With this notation the six (0,2) n, n' states can be written 
as products of functions in LR- and n-space. 



Rn' 
Rn 



= |tlR) |s„,„. 



(23) 



The tunneling Hamiltonian behaves as an identity op- 
erator in n-space, whereas in LR-space it acts as follows, 

^2p |^LR^^ ^ _y2t |t^^) , (24a) 

iJ^P |tLR) = -,/2t {\Tf) + JT^R)) , (24b) 

Hl^ \T^) ^ -V2t |T^R) , (24c) 

L^) = 0. (24d) 



iJ^P|S 



As a consequence of Eq.(24d), all the LR-antisymmetric 
solutions (the (1,1) states with a singlet JS'^'^)) are not 
affected by the inter-dot tunneling, i.e., they are blocked 
states. Combining Eqs. (15), (23) and (24a) it follows 



(21a) that the (0,2) solution 



Rn' 
Rn 



couples only to the even 



combination of the given (1,1) Slater determinants, 
1 



LR" 



) = 71 (|En') + |r:')) = |T^^) |Sn,n'> • (25) 



In the above mentioned space, the effect of the intcr-dot 
tunneling can be condensed into a 2 x 2 matrix 



rrovon ZP \_ ZP 



/2t 



Rn' 
Rn 



/2t 



and 



(26) 



LR" 



The maximum mixture between 

happens for £ = 0, where the avoided crossing occurs with 
a gap given by Aj = 2\f2t. Figure 3 shows the energies 



of the three two-particle states associated with the two 
single-particle states \n) and \n'). Six groups of 3 states 
with exactly the same energy gap appear, one for each of 
the possible n, n' non-equivalent pairs (see Fig. 1(c)). The 
four remaining states (that complete the 22 solutions) are 
the (1,1) LR- antisymmetric states that arise from Slater 
determinants constructed with the same single-particle 
state in both dots: 



IRti 
Ln 



) = |Si^«)(|n),|n),) 



(27) 



Using spin and valley singlet and triplets functions the 
states in Eq.(27) can be written (for the quantum dots 
of case A, B and C: ct is a good quantum number) as 
the four states with full polarization in spin and valley 

and 



cLRrnvalrpspm 



cLRrpvalrnspin 



spaces. I.e., 

In summary, there are ten LR-antisymmetric states 
that do not mix with the (0, 2) solutions. Six of them 
with n 7^ n' , as it is shown in Fig. 3, become quasi- 
degenerated with their associated n, n' LR-symmetric 
partners at the high detuning limit, i.e., when \£/t\ ^ 1. 
The six LR-symmctric states are connected by inter-dot 
tunneling to the corresponding n, n' (0,2) states as fol- 
lows, 



T^^) IS. 



^ 



|T^^) IS. 



(28) 



Obviously, as it is seen in detail in the next sections, the 
connectivity between (0,2) states and (1,1) states plays 
a key role in the return probability experiment. 



atoms in the dot and A{so the isotropic hyperfinc coupling 
constant. 

From Eq.(30) it follows that Eq.(29) with i = 3 
has zero coefRcients, as expected from absence of time- 
reversal symmetry. K i — Z term would correspond to 
a time-reversal symmetric spin-orbit interaction, which 
cannot originate from the hyperfine interaction. In the 
following sections we show that this apparent valley 
anisotropy of the hyperfine field have important conse- 
quences for the return probability experiment that man- 
ifest for some specific prepared states. 

The Hamiltonian of Eq.(29) and the variances in 
Eq.(30) are derived as follows. ^^ As in Ref. 25, we as- 
sume that the Fermi contact hyperfine interaction dom- 
inates. Furthermore, in order to minimize the number 
of special cases, we take the case with largest hyperfine- 
induced mixing, namely the case of isotropic coupling, ^^ 
even though there is some degree of anisotropy.^^ The 
hyperfine coupling Hamiltonian then reads 



Ht 



hf 



A;, 



•E 



CT ■ I/,s(5(r - R;, 



(31) 



where I is the unit cell index, s E {A, B} is the sublattice 
index and I;.s is the nuclear spin of the carbon atom at 
site His, being zero or a spin-1/2 operator for a ^^C or a 
^^C isotope, respectively. The normalized tight-binding 
eigenstates, characterized by the spin a and valley r in- 
dexes are 



(* 



^irKR,. 



T'^ll.. 



*I(Rm)xo 



(32) 



E. Hyperfine coupling 

Taking cr the spin operator of the confined electron in 
a graphene-based quantum dot, the hyperfine interaction 
has the following form:^^ 



7J;, = ^T,h«., 



(29) 



i=0 



The values of the h*^*' fields can be considered fixed dur- 
ing every single-shot measurement because the dynam- 
ics of the nuclei spins is much slower than the electron 
spin's. Within this approximation, the effect of the slow 
dynamics of the hyperfine field — on the average result of 
a experiment that is performed many times — is captured 
by noting that the hyperfine field components become 
random variables that follow zero-mean Gaussian distri- 
butions with variances^^ 



(0) 



,(3) 



= 2 
= 0, 



(i; 



= 2 h 



,(2; 



^ 2^L 



J ^ x,y,z, 



(30) 



where (Th = Af^^i'/{8NQo), with Nq^ being the number 
of atoms in the quantum dot, u the abundance of ^^C 



where the smoothly varying envelope functions ^^ (r) are 
eigenstates of the gapped Dirac equation in a confining 
potential. Taking matrix elements of these eigenstates 
with respect to the hyperfinc Hamiltonian then gives 
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*I(Rm)- (34) 



After performing averages over the nuclei fields, the 
Hamiltonian of Eq. (29) follows. The real (imaginary) 
part of Fi '~ having cosine (sine) factors generates the 
terms containing ti (T2). The two terms diagonal in val- 
ley Fl'J' , T ~ ±1 are equal and do not have oscillating 
factors, therefore, only operators proportional to tq sur- 
vive and the t^ component vanishes (leading to the im- 
portant property already mentioned). The factor of 1/2 

between the variances of the h^ ' and the h^ compo- 
nents comes from the averages over the cosine and sine 
squared. 

In a DQD, due to the different nuclei environments 
for each dot, it is convenient to work with the left /right 
homogeneous and the left/right inhomogeneous compo- 
nents of the hyperfinc field interaction (HFI) 
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In terms of these fields, the double dot single-particle HFI 
Haniiltonian can be written as 



H, 



ip 



E E /^feo + ^^'6 na„ (36) 



2—0 j^x.y,z 

which should be added to the Haniiltonian of the double 
dot in Eq.(7). 

1. Hyperfine field in single valley systems 

Here we briefly introduce the key features of the HFI 
found in single valley systems, e.g., GaAs double dots. 
The well-known relevant states and the action of the HFI 
Hamitonian in those systems is an important reference 
for comparison. Because the valley space is absent in 
GaAs, the single-particle Hamiltonian of the hyperfine 
field interaction is 
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(37) 



Due to the absence of the orbital degree of freedom 
there exist four (1,1) states, namely, the spin singlet and 
triplets. 
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The HFI Hamiltonian of Eq.(37), written in the latter 
basis, becomes 
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(39a) 



(39b) 

Sxy I— >■ v2 {Shx — \Shy) , xy h- > \/2 \hx — ihy) ■ (39c) 

One sees that only the inhomogeneous HFI is able to mix 
the LR-symmetric function (spin singlet) with the LR- 
antisymmetric (spin-triplets) functions. In particular the 
Shz component also conserves the total 5*2 and therefore 



oLRrpspm 



In 



is able to mix the ItJj-i^S'^p'") with the 

the next section we show that some specific situations in 
CNT double dots can be mapped to the Hamiltonian of 
Eq.(39a) or to a modified version of it. 



III. SIMULATING THE EXPERIMENT 

The return probability, P(ts), for a given evolving 
time, Tg, see Fig. 1(b), is obtained experimentally by av- 
eraging over a set (i = 1, 2, ..N) of single shot measure- 
ments with outcomes Oi^r^ , 



1 ^ 



(40) 



Each single-shot measurement consists of a gate-voltage 
cycle with five stages (see arrows and points "p" , "e" and 
"m" in Fig. 1(a) and the associated detuning values Sp, 
£c and Em in Fig. 1(c)): 

(i) Preparation, the DQD is prepared in the (0,2) region 
at point "p" . 

(ii) Separation, by applying a voltage pulse of length 
Ae and duration Ts, which is short on the scale of the 
HFI interaction {Ts ^ ^), but slow on the scale of the 
inverse tunneling energy, the initial (0,2) state is adiabat- 
ically moved to the point "e" deep into the (1,1) region. 
If the detuning is changed by Ae ^ t, the condition for 
adiabatic conversion is Ts » HAe/t"^. 

(iii) Evolution, the system is left to evolve at the point 
"e" during a time t^; in this stage the electron wave- 
function in each dot acquires a different dynamical phase 
due to the hyperfine coupling with the nuclei spins; in 
general, the system oscillates between the initial (1,1) 
wavefunction and other combination of (1,1) states. 

(iv) Joining, a voltage pulse brings the system back to 
the (0,2) region; with the same adiabatic condition as for 
the separation stage. 

(v) Measuring, at the "m" point a nearby charge sens- 
ing device determines the outcome Oi^r^ (1 or 0) of the 
single shot measurement; Oi,r^ is set to 1 only when the 
system has returned to an (0,2) configuration. 

Depending on the preparation protocol the initial (0,2) 
state can be different than the ground state. We cover all 
the possibilities by assuming a prepared state in an ar- 
bitrary superposition of the six possible zero-hyperfine 
(0,2) eigenstates, \^pp) = X^ti '^d^' (0'2)), where each 
|/, (0, 2)) state is one of the six Slater determinants given 
in Eq.(12). When averaging the outcomes Oi^r^ of a large 
number of single-shot measurements over an ensemble of 
random hyperfine field, the phases of a; average out and 
the resulting probability for a given initial state is simply 
given by 
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(41) 



where Pi{ts) is the return probability when starting in 
the (0,2) eigenstate \l, (0, 2)). Therefore, the return prob- 
ability of a general case can be evaluated by knowing the 
values I a; I and the return probabilities obtained when 
preparing the (0,2) zero-hyperfine eigenstates separately. 
For these reasons, in each of the DQD scenarios that we 
deal over the next sections, we focus on the behaviour of 
the six functions, Pi{ts). 

Due to its experimental importance we focus the dis- 
cussion on the saturation value of the return probability 
which is defined as 



P^= lim P{ts). 



(42) 



In addition, we define the decaying time, Td, as the time 
for which 



P{Td) = ^Poo + \, 



(43) 



see Fig. 1(b). 

In general, P{ts) may be different when the hyperfine 
field of the two dots follow Gaussian distributions with 
different rms values, a}^ and ct^. In fact, this is the case 
when the numbers of ^^C atoms are different for the left 
and right dots (see Eq.(30)). However, this would affect 
the transient of P{ts) but not the behavior for t^ — > oo. 
In this work we focus on the saturation values and on 
transient features intrinsic to DQDs in C-based systems. 
For this reason we choose aj^ = '^\ ^ <^h in the following. 



A. Numerical evaluation of the return probabilities 

We have developed a numerical simulation of the ex- 
perimental cycle outlined above. The time evolution 
inside the (1,1) region (point "e" in Fig. 1(a)) is gov- 
erned by the hyperfine field alone, i.e., it is assumed that 
t'^/Ae <^ at- Essentially, as discussed in the introduc- 
tion, we are working in the large detuning limit where the 
effect of the tunneling exchange is negligible and therefore 
we are able to obtain the lower bounds of the saturation 
return probabilities. For solving such a time evolution 



we start with the LR-symmetric (1,1) state. 
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even 



which is connected (see Eq.(28)) with the (0,2) prepared 



state 



Rn' 
Rn 



) ; this is then decomposed into the numerically 

determined eigenstates of the 16 by 16 Hamiltonian for 
the current hyperfine field realization (generated follow- 
ing the Gaussian distributions with variances given in 
Eq.(30)), and hence the time dependent (1,1) state can 
be computed as a superposition of time-evolved eigen- 
states. 

To calculate the return probability, we define an oper- 
ator that projects onto LR-symmetric states as 
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(44) 



The return probability after time t^ for a given hyperfine 
field realization r^ is then 



Pr^Ts) ^ (*r,(Ts)| ^LR-sym |*r.(Ts)) 



(45) 



The projection method assumes the joining stage (iv) is 
performed under the same adiabatic conditions as the 
separating stage (ii). After repeating this procedure for 
a large number of realizations, N^., the final return prob- 
ability P{ts) is obtained by averaging 
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(46) 



(i.e., no tunneling exchange), and that the return prob- 
ability after evolution can be computed by projection. 
In this high detuning limit — excluding for the moment 
the HFI effect — the LR-symmetric states can be consid- 
ered degenerated with their LR-antisymmetric partners 
(see Fig. 3 and Eq.(13)). Additional degeneracies are de- 
termined by presence of spin-orbit coupling and applied 
magnetic field. We assume that with finite spin-orbit 
coupling and/or applied magnetic field the hyperfine in- 
teraction only mixes states within the subset of quasi- 
degenerate states. 

Due to the LR-symmetry, since the evolution starts in 
a LR-even state, the degeneracies of the subspaces of evo- 
lution in (1,1) are never lower than two. In this paper 
degeneracies ne = 2,3,4, 8 and 16 appear. To determine 
the time evolution, we project the Hamiltonian to the 
(1,1) subspace connected with the prepared (0,2) state, 
the reduction of the system — if any — sometimes allows 
analytical treatment. Then, we focus on the form and 
the statistical properties (i.e., variances) of the surviving 
mixing terms of the HFI Hamiltonian. Those matrix ele- 
ments are obtained by elementary calculations; the (1,1) 
basis expanded using singlet and triplet functions in LR, 
valley and spin spaces (instead of the Slater determinant 
basis) is useful for a more direct physical interpretation 
of the hyperfine mixing terms. 



C. State counting estimation of Poo 

Here we introduce a scheme for estimating the re- 
turn probability based on a simple state counting ar- 
gument. The value obtained with the procedure that 
follows does not always coincide with the exact value. 
Poo 7 however it is useful for visualizing special features 
of the exact dynamics. In order to compute the state 
counting value, first, we find the number of degenerate 
(1,1) states, TT-e, connected with the chosen (0,2) prepared 
state under investigation. Second, we find the number of 
LR-symmetric states, ntQ2), connected with that (1,1) 
subspace of n^. states, Sg- Third, assuming a fully in- 
coherent mixing of the initial state with all states, for 
Ts — > oo every one of the He states should have a prob- 
ability l/rie of being occupied, therefore, the estimated 
return probability becomes 
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B. Analytical evaluation of the return probability 

For analytic evaluation of the return probability, we 
use the same set of conditions as for the numerical eval- 
uation, namely that the time evolution after the sep- 
aration stage is only governed by hyperfine interaction 



One can expect this to be a lower bound for Poo, be- 
cause under coherent evolution the system does not fully 
randomize and might therefore maintain a larger weight 
on the initial state which is connected. This is the case 
for example in GaAs double dots, whereas the estimation 
gives P^ = 1/4, the coherent evolution (averaged over 
many realizations) gives Poo = 1/3. 



IV. RESULTS 

1. Labeling of the energy levels 

With the exception of Sec. IV D (case D), in the fol- 
lowing, the single-particle eigenstates in each dot are 
|^,n) with |n) — |rcr), where a is the spin projection 
along the direction of the applied magnetic field taken, 
for convenience, along the ^-direction. The two-particle 
functions |S„_„') and |T„^„/) presented in Eqs. (22) can 
be further expanded in terms of tensor products in spin 
and valley spaces. The procedure is straightforward and 
states |S„.„') and |T„.„/) are found to be equivalent to 
a unique tensor product in valley and spin except for 
n, n' = K t, K' | and n,n' = K I, K' t in which case they 
are given by, 
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These states are particularly important in the presence 
of spin-orbit coupling. ^'^^■-'^^'^^ For the sake of readability 
in the following we sometimes omit the ket symbol, |. . . ), 
when referring to a product state of triplets and singlets 
in left/right, valley and spin spaces. 

As it is shown in detail in the next subsections, for 
cases A, B and C the sixteen (1,1) eigenstates for zero 
hyperfine are associated with 3, 9 and 10 energy levels, 
respectively, that sometimes cross each other or become 
degenerated. These situations change drastically the out- 
come of the return probability experiment. In order to 
present the correspondence between the (1,1) states and 
these energy levels, we divide the states into two classes. 
In the first class, as in Sec. II D, we take six groups, each 
one corresponding to a LR-even/LR-odd pair of states, 
i.e., |T^^S„,„') and \'S>^^Tn,n')- These six groups— that 
as shown in Table I, may belong to the same energy 



TABLE I. Energy levels and labels for the six LR-even/LR- 
odd partner states in the (1,1) configuration for cases A, B 
and C. We use the same labels for the connected energy levels 
in (0,2) to identify the prepared state. 
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TABLE II. Energy levels and labels for LR-odd full spin and 
valley polarized states in the (1,1) configuration for cases A, 
B and C. 
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level — are of great importance for the return probabil- 
ity experiment because they are connected to the (0,2) 
states, JT'l^Sn.n'), through the LR-even states. Apart 
from a global energy shift, these (1,1) energy levels and 
each associated (0,2) energy level have the same depen- 
dence with the parameters, therefore, we use the energy 
level labeling of Table I also for identifying the (0,2) pre- 
pared states. The second class of states consists of the 
four fully polarized spin and valley LR-odd states, i.e., 
\S^^) {\n)-^ 111)2)', they belong to the energy levels listed 
in Table II. These blocked states do not provide access 
to the (0,2) configuration neither at the preparation nor 
at the measurement stage, however, if their associated 
energy level crosses (or is degenerated with) an energy 
level with states of the first class, they become relevant 
for the effective hyperfine dynamics. 

The labeling introduced in Tables I and II for the en- 
ergy levels of case B and C — the less degenerated cases — 
is inspired by the following logic: "1" ,"c" and "h" denote 
low, central and high, respectively, referring to the valley 
characters of the states associated with the energy level. 
Similarly, the -I-, — and refer to the spin characters. 
Since two pairs of LR-partners states fall in level Cq, we 
distinguish them (in case C because Ago 7^ splits the Co 
level) by adding the subscript "H" or "L" for the high 
and the low energy LR-partner states, respectively. 



2. Labeling of the different physical situations 



In what follows the numbers Ur. and 



1(0,2; 



are useful 



for distinguishing different qualitative and quantitative 
dynamical situations. For this reason, we label each case 
using the two integers as no„,„ „, . As an example, with 

(0,2) 

the latter convention, the case in which the (1,1) sub- 
space is fully degenerated is to be labeled as "IGg" be- 
cause rio — 16 and n(o.2) = 6; in such a case six indepen- 
dent prepared states provide access to the full subspace 
of evolution. As another example, a crossing between the 
energy levels Cql and Iq (presented in Table I) implies a 
subspace of evolution with ric = 4 containing two LR- 
symmetric states (one associated with level Cql and the 
other with level Iq) that can be prepared, therefore, this 
situation is to be labeled as "42" . 

However, as we show below, in some cases two situa- 
tions with the same n^ and '^■(0,2) numbers but different 
HFI dynamics appear; in such cases we use a superscript 
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FIG. 4. Case A: Saturation return probabilities, energies and states for a valley degenerated DQD; Ago = Eorh = and the 
Zeeman splitting is Es—gsfJ^BB. The associated single dot single-particle spectrum is shown in Fig. 2(a). (Left panel) Eigenstates 
and spectrum of the (0,2) and (1,1) configurations; the degeneracy of each energy level at finite magnetic field is indicated. 
The six (0,2) states in levels v± and vq can be prepared; each one of them is connected by inter-dot tunneling, represented 
by double headed arrows, to a single LR-symmetric state in (1,1). (Right panel) Energy levels in the (1,1) configuration as a 
function of the magnetic field; as shown in the legend, P<x> is found to be 0.4, 0.50 or 3/8 depending on both the magnetic field 
value and the prepared (0,2) state. 



"s" or "n" to distinguish them. The "n" superscript is 
reserved for cases in which the effective HFI Hamilto- 
nian has a zero matrix element between any of the LR- 
symmetric states in the evolution subspace, Tg^S„^„'), 
and its LR-antisymmetric partner state, S'^^Tn.n'). 



A. Case A: Aso = and _Borb = 

The single-particle spectrum of a single dot (in the ab- 
sence of HFI) is shown in Fig. 2 (a) as a function of the 
magnetic field for this case. The magnetic field, Bz, can 
point in any direction as long as diamagnetic effects are 
avoided and, therefore, the valley degeneracy is not lifted; 
for nanotubes this impose that By = 0, whereas for a 
graphene quantum dot the magnetic field must be in the 
plane. ^^ The energy levels for the (0,2) configuration are 
only three, v^^ with s^ = ±, (as listed in Table I), with 

energies E^^'^ — SzEg] this follows from Eq. (12). In the 
left panel of Fig. 4, we show the states and label them in 
the singlet/triplet notation for spin, valley and left/right 
space. Also the degeneracies at finite magnetic field are 
given for each group. The connections of (0,2) and (1,1) 
states by tunneling are shown with double headed arrows, 
such that only left/right symmetric (1,1) states are con- 
nected to their (0,2) partners with same spin and valley 
quantum numbers. 

The right panel of Fig. 4 shows the energies of the (1,1) 
states as a function of the magnetic field. The saturation 
return probability is also given in the figure as tags to 
the energy versus field lines. These symbols mean that 
if the (0,2) system is prepared in vq or y± the return 
probability is 1/2 and 3/8, respectively, while for B = 0, 
the return probability is 0.40. Below we explain each of 
these cases in more detail. 

The Case IGg (in Fig. 4) is for zero magnetic field, 
leading to the largest possible degeneracy, n^ = 16. All 



the (0,2) eigenstates are connected to this subspace and 
hence the state counting estimation for the saturation re- 
turn probability is P^'^^*^ =6/16 = 0.375, which, however, 
does not coincide with the actual result P^*^ =0.40. The 
result is found irrespective of which of the six (0,2) states 
is prepared. The increment of the return probability, as 
compared to the state counting estimation, is similar to 
the GaAs situation, where the state counting gives 1/4, 
whereas coherent evolution gives return probability 1/3. 
Interestingly, in both cases one goes from P^ to the cor- 
rect result by subtracting 1 from the denominator, since 
6/(16-l)=0.4. 

We have also studied the final state resolved satura- 
tion return probabilities and observed that it is not the 
same for the six (0,2) states. It is more probable to re- 
turn to the (0,2) state that was prepared, and this is 
more pronounced when the prepared state is one of the 
two nonzero-spin states (the states t^^S™'T^'" ). In 
Fig. 5 (a) we show that the decaying time r^ (see Fig. 1(b)) 
here depends on the prepared (0,2) state, being longer 
for the nonzero-spin prepared states, 0.185h/(Jh, while 
for the remaining four prepared states it is 0.1A9h/ah. 
We discuss below, when introducing cases 8| and 4", the 
reason of such an asymmetry. 

It is worth to introduce here the structure of the hy- 
perfine field interaction Hamiltonian, Hf^^^^, using the 
sixteen (1,1) basis states presented in Fig. 4 and Ta- 
bles I and II. In the octahedral representation given in 
Fig. 5(c) each vertex stands for a LR-symmetric state and 
its LR-antisymmetric partner state. Therefore, the ini- 
tial (1,1) state at Ts = has weight in one of the LR- 
even states located on those six vertices. The links in 
solid lines between vertices and the links in dotted lines 
with the remaining LR-antisymmetric states represent 
nonzero matrix elements of the HFI. Each link stands 
for both homogeneous and inhomogeneous elements of 
the hyperfine field. The absence of matrix elements be- 
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FIG. 5. Scheme of the complete hyperfine Hamiltonian (right panel) and the return probabilities P{ts) for case A (left panels), 
i.e., Aso = £^orb = 0; see Fig. 4. (a) Zero magnetic field, type IGe dynamics, all prepared states give Poo =0.40 but P(ts) decays 
faster for total spin zero prepared states, (b) Finite magnetic field. Type 84 dynamics is found when a zero spin state is 
prepared, Poo =0.50. For nonzero spin prepared states one has type 4i dynamics and P{rs) follows Eq.(54) giving Poo =3/8 
(this situation also appears for the DQDs with Ago 7^0 of case C, see Fig. 8). (c) Scheme of the full 16 by 16 HFI Hamiltonian, 
HI^''\ we label the states according to the ten energy levels they belong in case C. The six octahedron's vertices represent 
LR-odd/LR-even partner states of Table I; note that the HFI matrix element between these partners, m, is zero for the c_|. and 
the C- vertices. Outside the octahedron the four LR-odd states of Table II are represented (for simplicity, only two are shown). 
The links between nodes represent the existence of at least one nonzero matrix element between the states. Depending on the 
parameters — for cases A, B and C — the ten energy levels (c±, cqh, cql, ho, lo, h± and 1±) can move together and/or cross each 
other, leading to the eight different effective dynamics involving less than sixteen (1,1) states. 



tween states in opposite vertices of the octahedron follows 
from the trivial selection rule that forbids two single- 
particle quantum numbers from being changed simulta- 
neously; note that the two-particle HFI Hamiltonian is 
Hf^i.K = li ® (^/.!'lr)2 + (^h!LR)i ® I2. This rule also 
justifies that each one of the four IS'"^) {\n)^ W)-^) states 
(in the figure nodes h± and \± represented outside the 
octahedron, for simplicity only h_ and 1+ are shown) 
mixes only with the three LR-even/LR-odd partners — 
rules given in the figure — that have one of the two elec- 
tron in the \n) state. 

In Fig. 5(c) we have used the labeling of the energy lev- 
els given in Table I for case C (see Fig. 8), c±, Cqh, Col, 
ho, lo, h± and 1±. We show below that — depending on 
the parameters in cases A, B and C, as some of these ten 
energy levels move together and/or cross each other — 
there are eight different effective HFI dynamics involv- 
ing restricted subspaces of evolution. For example, in 
the presence of spin-orbit coupling (case C below) the 
IGg case becomes irrelevant because the hyperfine field 
(characterized by an energy scale, cr^, which is smaller 
than Ago) is unable to mix all the sixteen states and 
therefore the physical situation is better captured by an- 
alyzing smaller subspaces of evolution depending on the 
prepared state. 

The Case 8|. A finite Zeeman energy is applied much 
larger than ah- The initialized (0,2) state has zero spin 
and therefore belongs to the energy level vg with degener- 
acy ?T-(o,2) =4. The corresponding (1,1) subspace in which 
the system evolves has a degeneracy of no = 8, with four 



LR-symmetric states. In this case, the state counting 
estimation for the saturation return probability, being 

Poo' * = 4/8 = 0.5, coincides with the calculated exact 
result. Moreover, as shown in Fig. 5(b), the values of Pqo 
and the shape of P{ts) do not depend on which of the 
four (0,2) states is being prepared. Finally, we find that 
the decaying time, here Td sa 0.238/i/crh, is larger than 
both possible decay times for the zero-field case, which is 
expected because for the 165 case the HFI has the ability 
to mix the prepared state with eight extra states. 

The independence on the prepared state can be un- 
derstood from the matrix elements of the Hamiltonian 
in the reduced Hilbert space by using the singlet/triplet 

functions. The HFI term proportional to 5hz mixes one 
by one the four zero-spin pairs of LR-symmetric and LR- 
antisymmetric states as follows, 

(S^«TKt,K';| i^^^LR It^^'Skt^k';) = 2<5/i("\ (49a) 
'/.,LRrrSKi,K't)=-2<5/if,(49b) 



/nLRrp I rr2p 

\0 i-'Ki,K'^\-tl^. 



cLRrjivalrpspin 







cLRrjivalrpspin 







^^XR |T|^«T™'S^P'") = 2(5/li"\ (49c) 

^ft^LR JTjj'^T^'S'P'") = 25h^°\ (49d) 



The fact that the subspace S'o* is composed solely by LR- 
symmetric states and their LR-antisymmetric partners is 
a key difference to the other cases where the state count- 
ing argument fails, because in those cases more states are 
available. 

An illustration of all the matrix elements Hc^ can be 
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extracted from the full Hamiltonian representation in 
Fig. 5(c), by considering solely the four octahedron ver- 
tices in the plane (because the vg level contains all the 
states in nodes Ig, hg, Cqh and Cql) and excluding the 
c+ and c_ vertices. In accordance with Eq.(49), the 
HFI matrix element between the included LR-partners, 
m in the figure, is nonzero. The remaining matrix ele- 
ments (links in the figure) connect each LR-synimetric 
state with other two LR-symmetric states and their LR- 
antisymmetric partners. All those matrix elements follow 
Gaussian distributions with the same variance and thus, 
the average result is independent of which Tg^ state is 
initialized. This symmetry is also responsible for the dis- 
tribution of the total saturation probability (1/2) evenly 
among the four Tg^ states. 

The case 4" (Fig. 4) is also for E^ ^ ah but here the 
initial state belongs either to the energy level v+ or the 
v_. In Sec. IV C the same situation is found for a par- 
ticular value of the parallel magnetic field in a nanotube 
with Ago 7^ 0. Here we present an analytical derivation 
of the return probability in the reduced Hilbert space, 
which we have also confirmed by a full numerical evalu- 
ation. We find that the saturation return probability is 

p4" =3/8 = 0.375. 

We restrict the analysis to the states of the energy 
level v+, and equivalent results for level v_ follow by 
symmetry. The system is prepared in the t!i^S™'T^'" 
(0,2) state, and after separation the initial (1,1) state is 
therefore T^i^S^'r/" 
4 and the subspace is composed by 



kj -L+ • The degeneracy of v+ in (1,1) is 



|1,45') = tIj-'^S^'T^p'", 12,45') - s^^T^'r/", 
13,45') = s^Rt^'T^p'", |4,45') = s^^T^'^/". (50) 

There is only one LR-symmetric state (i.e., n(n2) = 1) 
and it is a spin-polarized valley singlet. By restricting 
the full HFI Hamiltonian to the latter subspace we get 
the eff'ective 4x4 Hamiltonian: 



ffo° ^2U?H.+H?^^\ 



(51) 



where I4 is the 4-dimensions identity matrix and H^'^^^ 
is the 4x4 Hamiltonian presented of Eq.(39a) after the 
following replacements, 

Sz^Q , ZH^O , (52) 

Sxy ^ V2 f(5/iW _ iSh'-^A , xy^V2 f /i^ _ [hi^)\ . 

Note that the matrix elements between the nonzero spin 
LR-symmetric/LR-antisymmetric pairs is 



Eq.(53) is due to the absence of terms proportional to 
the T3 operator in the HFI Halmiltonian of Eq.(36); as 
discussed above, such terms would preserve time reversal 
symmetry and thus cannot originate from the hyperfine 
field. 

It is instructive to compare He ' with the effective 
Hamiltonian for zero field in GaAs double dots (see 
Sec. HE). The null matrix elements are those that in 
a 2-dimensional electron gas DQD would translate to the 
Shz and hz components. In our system, the effective 
hyperfine field operates on valley-space (all members of 
the subspace in Eq. (50) share the same two-electron spin 
function) and there is no mixing of the valley singlet, 
S™', with the valley triplet, T™', nor any splitting of 
the T™' states. Hence, such a valley double dot maps 
to the typical "spin only" double dot but with an ef- 
fective 2-dimensional hyperfine field instead of the usual 
3-dimensional HFI. 

In Appendix A we give the details of a standard analyt- 
ical procedure for solving the time dependence dictated 
by He ' and averaging the result over the Gaussian fields. 
We obtain the return probability, shown in Fig. 5(b), 

P^"(ts) = C2B {rs<Jh/h) C2B {r^af/h) 

+252D (r.^rL/fi) 52D {r^af/h) , (54a) 

-2^'crfi {y2x\ , (54b) 

C2D {x)^l- 2S2D (x) , (54c) 



^^^("^^^yr^ 



where the imaginary error function is erfi(x) =erf (ia;)/i = 
^~^ J- dyexp{y^). As already mentioned, the satura- 
tion return probability is 



lim S'2D(a;) 



1 



^ 00 



3/8 = 0.375 



(55) 



-,sc,4? 



For the present case, the state counting value is Po< 
1/4, and therefore this is another example of a case in 
which Poo is greater than the state counting value. Fur- 
thermore, here the substraction of 1 to the denominator 
in the state counting estimation docs not provide a cor- 
rect result as it does for the zero magnetic field cases in 
a GaAs DQD or type IGg dynamics. Finally, applying 
Eq.(43) to the return probability P^i (r^), we get the de- 
caying time Td ~ 0.3218h/ ah, which is the slowest time 
we find for case A; i.e., in C-based double dots with un- 
broken valley degeneracy. 



B. 



Case B: Aso = and £orb/0 



oLRrpvalrpspin 



Jj2p 
-"/i,LR 



rpLRgvalr[.spm ^ ^ q 



(53) 



This expression, together with the four zero spin cases 
presented in Eq.(49), completes the effective mixing be- 
tween the six Tg^/S'"^ partners; Fig. 5(c) shows the six 
pair of states together emphasizing that the mixing, m, 
is zero for the two pairs presented here. Remarkably, 



In Fig. 2(b) we show the single-particle spectrum of a 
single dot as a function of the parallel component of the 
magnetic field, P|| . Assuming that the total magnetic 
field is along the z-direction at an angle (j) from the tube 
axis, we get Py = B^ cos <j>. The valley degeneracy is 
also lifted in contrast to the case without orbital mag- 
netism studied in Sec. IV A. In the left panel of Fig. 6 we 



13 



(0,2) configuration 

6 states 





ho 1 






-e 






> 


c+ 1 

Co 2 

' c. 1 


•e 






\ 


,lo 1 



ry LRry ValoSpill 



■- KT.K'i 




Ki.K'T 
\T-' LR c valT-- spin 

^T-- LR T-- valo spin 



Interdot 
tunneling 

■* ► 



(1,1) configuration 

LR symmetric h j LR-anti symmetric states (10) 
states (6) 



K 



'yLK'yvalQS] 
'-T->LRQvaI'-T->s| 



'o2 
h 1 



2£. 



oLR'yval'yS 



To^^S 



Kt.K'i 



rpLRc 
T-'LRcvalT-' 



T-'LRT-'valQS] 



KT.K'i 




Ki.K't 



^ oLRT^valT-'S 

oLRT-'valT-'spin 

(1,1) States with 

the same single 

particle state in 

the two dots: 



1 



2E, 



Q LRt-- valT-- spin 



oLRT-'vaJT-'Si 



oLRT-'valq-'Si 




FIG. 6. Case B: Saturation return probabilities, energies and states for a nanotube-based DQD without spin-orbit coupling 
in a magnetic field with a parallel component (in the figure we take -Borb > Eg), i.e., Aso = 0, -EorbT^O and Eb^O; the single 
dot single-particle spectrum is shown in Fig. 2(b). (Left panel) Scheme of the (0,2) and (1,1) energy levels and the associated 
states' singlet/triplet characteristics — in spin, valley and left-right space; the level degeneracies for finite magnetic field are 
indicated. Each one of the six (0,2) states is connected by inter-dot tunneling to a single associated (1,1) state. For the (1,1) 
charge state the six LR-symmetric and the ten LR-antisymmetric states are distinguished. (Right panel) Energy levels in the 
(1,1) configuration as a function of Bf, dashed levels contain only blocked states and therefore cannot be prepared. As shown 
in the legend. Poo is found to be 0.40, 1/2 or 1 depending on both the magnetic field and the prepared (0,2) state. 



present the spectrum for the (0,2) and (1,1) configura- 
tions including the level degeneracies at finite magnetic 
field and the corresponding states. In the figure, we take 
a ratio ii^orb/^s — C^gonh/ Qs) cos (j) bigger than one; which 
is the typical case for fully parallel magnetic field. ^"^ The 
five energy levels of the (0,2) configuration-c± , Cq, Iq and 
hg-follow from Eq.(23) and are given in Tables 1 and II. 

For zero field all the degeneracies remain and the situa- 
tion reduces to the IGg case already analyzed in Sec. IV A. 
At finite B\\ , the return probability experiment can have 
two different behaviors as indicated in the right panel 
of Fig. 6. When one of the four spin-zero states is pre- 
pared, we get Poo = 1/2, whereas a nonzero spin prepared 
state is unaffected by the hyperfine interaction and thus 
P^i (ts) = 1. We explain these two cases in the next para- 
graphs. 

In the case labeled 2\ the initial state belongs to 
the ho, lo or to the Co energy levels. The degeneracy 
for the (1,1) states in Iq or in ho is 2, and n(o.2) ~ 1, in 
each case corresponding to a pair of LR symmetric/anti- 
symmetric states. On the other hand the Co level has a 
4-fold degeneracy and two LR-even states are connected 
to (0,2); however, the structure of the hyperfine interac- 
tion Hamiltonian allows us to treat the subsets of double 
degenerated states, Col and Cqh presented in Table I, as 
two independent pairs with rio = 2 and n(o,2) = 1- This 
follows from the selection rule introduced above when 
explaining the absence of matrix elements for states at 
opposite vertices of the octahedron in Fig. 5(c). 

Here the state counting prediction for the saturation 

sc '2^ 

return probability (being Poo' ^ — 1/2) coincides with 
the exact result. The 2x2 effective Hamiltonian that 
describes the four cases follows from Eq.(49) and it is 



given by 



H, 



2'^ 



26h 

,(0) 



(0) 



(56) 



^ 26h>i" 
In Appendix A we derive that the return probability is 



P'^{rs) = 



(57) 



which is shown in Fig. 7. From Eq.(43) one readily ob- 
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FIG. 7. Return probabilities P{ts) for case B, i.e., zero spin- 
orbit coupling but coexistence of Zeeman and orbital effects 
due to the parallel component of the magnetic field (see labels 
in Fig. 6). Depending on the prepared state the return prob- 
ability have two behaviors: (i) type 2° — with Poo = 1 — when 
one of the two nonzero-spin states are prepared or (ii) type 
4" — with Poo ~ 1/2, see Eq.(57) — when one of the remaining 
four states is prepared. 



14 



f^(0,2) configuration 



(six states) 



1- Tl-'^S 



< 



Ki.K't 



rrn LRrjH valQspin 






LRcvalrT-ispiti ^ 



lo 1 



1 rT-iLRQValrT-i; 

V 

1 rT-iLRc valrT-ispin 



^rT-1 LRr-p vale spin 



T'-'^S. 



Interdot 
tunneling 



(1,1) configuration (sixteen states) 



LR symmetric 
states (six) 



h 1 



ryiLRc _ 



ry LR ry valo Spin 



ri-i LRq valrT-1 spin 
r-p LRq valr-p spin 



ri-i LR rT-1 valo spin 



ri-iLRc 
^0 ^Kt,lC'i 



OH 2 

1+ 1 



LR-anti symmetric states (ten) 

oLR'y valrT-ispin ■ 

,S-T, ■ ^' ^ 



Ki.Kt 



2(/?„, 



LRrjHvalrTHSpin 




oLRrTH valrjH spin 

Four states 
with iden- 
tical single 
particle states 
occupied in 
the 2 dots. 



LRrjHvalrTHSpin 



; cLR^val^spin 



Pec - 84",'H'l/2 ; ,n' I , 4n^^j/2 ' 4° ,^3/8; 




FIG. 8. Case C: Saturation return probabilities, energies and states for a nanotube-based DQD with spin-orbit coupling in 
a parallel magnetic field, i.e.. Ago 7^ 0, -Eorb 7^ and i?s 7^ 0; the single dot single-particle spectrum is shown in Fig. 2(c) (Left 
panel) Scheme of the (0,2) and (1,1) energy levels and the associated states; the level degeneracies for finite magnetic field 
are indicated. Each one of the six (0,2) states is connected by inter-dot tunneling to a single associated (1,1) state. For the 
(1,1) charge state the LR-symmetric and LR-antisymmetric states are distinguished. (Right panel) Energy levels in the (1,1) 
configuration as a function of the magnetic field; dashed levels contain only blocked states and therefore cannot be prepared. 
As shown in the legend — for the fifteen level crossings at Bi, with i = 0, . . . ,8 and the six generic situations at any non-crossing 
field, Bnc — the return probability presents seven different types of dynamics and Poo is found to be 1/3, 3/8, 1/2 or 1 depending 
on both the magnetic field value and the prepared (0,2) state. 



This type of dynamics is equivalent to the situa- 
tion for a GaAs double dot in large Zeeman field, in 
which case only the spin triplet S^^Tq'''" and the spin 
singlet Xq^S'^p™ are — neglecting tunneling exchange — 
degenerated; in that system the z-component of the hy- 
perfine field (see Eq.(39c)) is responsible for the mixing 
of these two states. ^"'^'"'^^'^'' In the present case an equiv- 
alent physical situation is realized for four out of the six 
possible prepared (0,2) states. 

In the case 2" the prepared state belongs to one of 
the two (0,2) energy levels c_|_ and c_, where the states 
are the spin-polarized valley singlets: 



rpLRgvalrpspin j^^ j^^^j ^^^ ^ ^j^j^ ^^ ^ ^ 



(58) 



After separation the system is initialized at the associated 
LR-even state, T^j'^S^^^'T^f", which is degenerated with 
its LR-odd partner, sLR-T^^T^p". As shown previously 
(see Eq.(53)), the hyperfine interaction does not have 
matrix elements between these two states. The return 
probability is therefore Tg -independent and equal to one 



P^Hr,) = l- 



(59) 



This somewhat counterintuitive result has been reported 
previously in Ref. 24, where the same saturation value is 
found in silicon double dots — taking into account a val- 
ley conserving HFI — for some particular prepared states 
that, as here, the hyperfine field is unable to mix with 
other states. 



C. Case C: Aso/0 and Sorb/ 

Here we consider a double dot based on the single dot 
presented in Fig. 2(c). Figure 8 shows the energy levels 
and states of the (0,2) and (1,1) configurations, which 
follow from Eqs.(12) and (13). The right panel depicts 
the ten levels of the (1,1) spectrum as a function of the 
total magnetic field which is parallel to the tube axis. In 
the figure, we use solid lines for the six energy levels hav- 
ing LR-even states, those states are accessible from (0,2). 
By preparing a state in those levels the return probabil- 
ity experiment can be performed; at the evolution stage, 
it becomes accesible to the HFI each and every state be- 
longing to others energy levels if they are degenerated 
with the energy level that holds the prepared state and, 
therefore, level crossings must be studied. As it is shown 
in Fig. 8, it is intrinsic to this case the existence of level 
crossings for zero magnetic field and also at finite values 
of the magnetic field. Finite magnetic field crossings (in- 
volving at least one energy level with a LR-even state) 
occur at 

Bi = 



B;: 



B., 



AiB(4gorb+gs)' ""^ /^B(4gorb -gs) ' 




Aso g _ A,o 
AiB(2gorb + 2gs)' ^ 2/iBgorb ' 


(60b) 


Aso g Aso 


(60c) 


MB(2gorb + gs) ' ° AiB(2gorb " 2gs) 
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B7 = 



MB(2gorb -gs)' 



B8 = 



A^Bgs 



(60d) 



these being positive values of the magnetic field given 
that we have assumed i?orb > ^-s and therefore 2(7orb > 
(7s- We identify twenty-one different situations, the six 
non-crossing cases and the fifteen level crossings. In the 
right panel of Fig. 8 we label each case according to its 
type of dynamics and we give the value of Poo in the 
legend. We describe all those situations below. 



1. Zero magnetic field 

As shown in Fig. 2(c) in the single-particle descrip- 
tion, the spin-orbit coupling breaks the 4-fold degener- 
acy resulting in two Kramers doublets; the lowest energy 
Kramers doublet in the quantum dot ^ (L or R) con- 
sists of the pair of time- reversal states |^K t) and |^K' ^), 
whereas the highest energy Kramers doublet groups the 
states I^K' t) and ICKJ,). With double occupation of 
the right dot — states and levels shown in Fig. 8 — for the 
state at level Cql (cqh) the two electrons occupy the low- 
est (highest) energy Kramers doublet in dot configuring 
the non-degenerated ground (highest excited) state of the 
(0,2) configuration. Right at the middle energy between 
the last two states the levels c+ , c_ , Iq and ho are degen- 
erated: they correspond to the four (0,2) states with one 
electron in each of the Kramers doublets. 

For the type 4| presented in Fig. 8 the return proba- 
bility behaves identically when the prepared state is the 
(0,2) ground state (level Cql identified below by cr = -|-) or 
the (0,2) highest exited state (level Cqh identified below 
by (J — —). After separation the LR-even state becomes 
a member of an (1,1) evolution subspace with nc = 4, the 
four states are 

llj^^i"') ^ Tq Skct.k'ct , |2,4^,(t)=S ■Tk<j,k'ct , 

|3, A\,<7) = sLRt™'T^/", |4, A\,cf) = sLRx^ix^pi", (61) 

where in the spin triplets the subindex a is to be in- 
terpreted as + or — instead of ^ or ^, respectively. In 
this space the effective Hamiltonian i/c,cr is analogous to 
a double dot in GaAs described by the Hamiltonian in 
Eq.(51). The equivalence with the components of the 
hyperfine field given in Eq.(38) is as follows 

SK ^ aSh'-^^ , h, ^ ah^°^ , (62a) 

Sh^ ^ 6hi^'> - cr6h^^^ , h^ h^ hi^'> - ah^^^ , (62b) 
Shy ^ Sh'-/^ + aShl^^ , hy ^ /i(,2) ^ ^^a)^ (g2c) 

From Eq. (30) we obtain that the components of the ef- 
fective hyperfine field 

hj — hj + Shj , h^ ~ hj — 6hj for j — x, y, z , (63) 

follow Gaussian distributions with identical standard de- 
viations (748 =\/2ah- Therefore, the situation is mapped 



exactly to a spin-only double dot at zero field and the 
return probability is then 

P^Hr.) = ^ (1 + 5 {rs<yii/h)+g'' {T^ai^Jh)) , (64) 

where g{x) = e^^^ (l — 4a;^). A derivation of this re- 
sult is presented in Appendix A. This well-known shape, 
shown in Fig. 9(a), leads to Poo — 1/3 and so the state 
counting estimation, which is this case is 1/4, fails. Fi- 
nally, for this shape the decaying time defined in Eq. (43) 
is Td^O.imihh/ah- 

In the case 8" the prepared (0,2) state belong to the 
energy level c+, c_, ho or lo. At Bo = the four levels 
have the same energy and therefore, after the separation 
stage, the evolution subspace in (1,1) includes their four 
LR-even partners and the associated LR-odd states, i.e., 
77,0 = 8. For all these prepared states we find Po^ =0.50, 
but the shape of P^'^{ts), shown in Fig. 9(a), depends on 
the prepared state. The decaying time is Td«0.301/i/CTh 
when the prepared state belongs to the level c+ or c_ and 
Td K, 0.201fi/(T/j when the prepared state belongs to the 
level ho or Iq. From the full HFI Hamiltonian in Fig. 5(c) 
one can visualize the Hamiltonian of the hyperfine field 
for this subspace, ife * , by considering states at the oc- 
tahedron vertices excluding Col and Coh vertices. The 
difference in the decaying times arises because, as shown 
in Eq.(53), there are no direct matrix elements of the HFI 
Hamiltonian between the nonzero-spin LR-even/LR-odd 
partners; i.e., to = for c_|_ and c_ in Fig. 5(c). 



2. Finite magnetic field 

There is a set of situations — including any non-crossing 
value of the magnetic field, B^c — in which the return 
probability is analogous to the cases without spin-orbit 
coupling introduced in Sec. IV B. This happens for the 
following situations (see crossings at Fig. 8): 

(i) When the prepared state belongs to the energy level 
lo at Bi (or ho at B2) because T™' and T™' are not 
mixed by the HFI. The return probability is type 

2\ (Poo = 1/2). 

(ii) When the prepared state belongs to the energy level 
c_ at B3 (or at Bg) because T*!^'" and T!f"" are not 
mixed by the HFI. The return probability is type 

2? (Poo = l). 

(iii) When the prepared state belongs either to the en- 
ergy level CoH (one electron in \K \) and the other 
in \K' t)) at B5 or, to level Col (one electron in 
\K^) and the other in \K' i)) at B7. At B5, the 
level Coh crosses h+ (both electrons at the \K f) 
state) and, at B7, the level col crosses 1+ (both 
electros at the \K' f) state). The hyperfine field 
does not introduce mixing at these two crossings 
because in each one of them the LR-odd state have 
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one gets Poo = 1/3 (type 4f ) or Poo = 1/2 (type 84). In the 
latter case the probability decays slower if the prepared state 
is a spin polarized valley singlet, (b) All cases give Poo = 1/2. 
For _Bz=B4 (type 4|) at the crossing (4,1) (or at the crossing 
(4,111)) P{t3) is independent of the prepared state. On the 
other hand (see text regarding the type 4^ case) at B^ — Bs 
the shape of P{rs) differs depending of which state in the 
crossing (8,1) (or crossing (8,11)) is prepared. At the crossing 
(8,111) (the LR-even state at ho is prepared) the behavior is 
the type 4| presented in panel (a). 

two single-particle quantum numbers different than 
the LR-even/LR-odd pair of states (see scheme of 
the full HFI Hamiltonian at Fig. 5(c)). The return 
probability is type 2f {Poo — 1/2). 

In addition, there are two special situations that lead 
to Hamiltonians already presented (see right panel in 

Fig.8): 

(a) The crossing (4,11) is relevant for the prepared 
state belonging to the energy level c+ at B4, which 
crosses the levels h+ and 1+ simultaneously. The 
four states in the evolving subspace have spin T^"" 
function, and the system works as a valley double 
dot with type 4" dynamics (see the Hamiltonian 

He^ in Eq.(51)), therefore, the return probability 
follows Eq.(54); i.e.. Poo =3/8. 

(b) The crossing (8,111) is relevant for the prepared 
state belonging to the energy level ho at Bg, which 
crosses h_|- and h_ . The four states in the evolving 
subspace have a valley T™' function. Since only 
terms in the HFI proportional to the tq operator 
can mix them the situation can be mapped to a 
non- valley degenerated (i.e., spin-only) DQD as in 
GaAs. This is achieved by the following replace- 
ments in Eq.(39c) 



This is a type 4| behavior and, as in crossings (0,1) 
and (0,111), the effective components follow a Gaus- 
sian distribution with standard deviation v2cr/i. 
The return probability is given in Eq.(64) and the 
saturation value is Pqo = 1/3. 

Finally, two double degenerated levels cross each other 
in the four remaining situations; namely, cases labeled as 
4| at (4,1) and (4,111), and cases labeled as 4| at (8,1) and 
(8,11). In each crossing the subspace of evolution has four 
states and two of them are connected to (0,2) and there- 
fore any of the two can be prepared. We find Poo =0.50 
for all the cases. However, as shown in Fig. 9(b), the 
dynamics in these two classes of crossings (see the right 
panel of Fig.8) are different: for type 4| the decaying time 
is, Td«0.2585/i/cr/i, independently of the prepared state, 
while, on the other hand, for type 42 the decaying time 
is Td « 0. 2327/1/ (T/i if the prepared state belongs to the 
level Cqh or Cql, or, Td^0A366h/afi if the prepared state 
belongs to the level c_ or C-).. The independence on the 
prepared state found for type 4| is justified by noting that 
its dynamics its governed by the symmetric Hamiltonian 
for type 8| — introduced in case A — in a smaller subspace 
that preserves its original symmetry. On the other hand, 
for type 43, the strong dependency of the transient on 
the prepared state arises because the crossings (8,1) and 
(8,11) involve the LR-even spin polarized valley singlet 
states at level c+ or level c_ . Those states dephase more 
slowly because, in accordance with Eq.(53), the HFI ma- 
trix elements with their LR-odd partners are zero. 



D. 



Case D: Aso/0 and Bj_/0 



Here, in contrast to cases A, B and C, the spin pro- 
jection along the direction of the magnetic field is not a 
good quantum number. The perpendicular field, B±, in- 
troduces a Zeeman energy Eg — gfiBB± and zero diamag- 
netic effects. Due to the competition between the Zee- 
man interaction and the spin-orbit coupling the single- 
particle and single dot problem has eigenstates with spin 
projection in the plane generated by the tube axis and 
the direction of the magnetic field. In the following the 
tube axis is chosen along the z-direction and the mag- 
netic field is applied along the a;-direction. The solutions 
are 



|K,+; 
|K',+; 
|K,-; 

IK',-' 



= |K)®(cos||i)+sin||t)), (66a) 

= |K')®(cos||t)+sin||;)), (66b) 

= |K)®(cos||t)-sin||;)), (66c) 

= |K')®(cos||;)-sin||t)), (66d) 



where 77 



arctanE's/Ag 



Shj i-> Sh 



(0) 



hj 1—^ h 



(0) 



for j = x,y,z. 



The eigenenergies E± — 
-2 V "so I P^s ^^^ shown in Fig. 2(d). From hereon we 
use the doublet index, d = ±, to identify the two dou- 
(65) blets. The tunneling Hamiltonian is still diagonal in this 
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basis, which means that each (0,2) state mixes with only 
one (1,1) state and the tunnehng energy gap is A^ = 2-\/2i, 
as before. 

We now show that the perpendicular field situation re- 
duces to a modified version of cases already considered in 
the paper. When the single-particle states are well sep- 
arated on the scale of the HFI only the matrix elements 
of the HFI Hamiltonian between the d = -I- solutions or 
in between the d = — solutions enter. The 2x2 effective 
hyperfine field for the dot ^ and the doublet d can be 
writing as 



m^' = 



e.d_d 



hi-''ai''' + hl'^al-'' + hf''ai'^), (67) 

where £, — L, R and cr^' are Pauli and identity matrices 
in the doublet space. The coefficients are given by 



hi''' ^ h^^^dsmf^, 


(68a) 


hi'^h^^l + h^^ldcosv, 


(68b) 


'^^'' = <J-'^2dcos77, 


(68c) 


hi''^ = -h^^Ucosr,. 


(68d) 



The Hamiltonian Hf^' is equivalent to a spin in a Zee- 
man field, plus an energy shift, Hq , which is irrelevant 
for the dynamics of the 2-level system. The values of the 
effective field components depend on the angle 77, i.e., on 
the external perpendicular magnetic field. Using the vari- 
ances of the hyperfine components h^^l (for j — x,y,z and 
1 = 0, 1,2) given in Eq.(30), it follows that the variances 
of the effective components in Hf' arc: 



H^' 



i^r 



aid 



?f^ (/^l'i)n = 2a^os^. 



cos^ rj) , (69a) 
(69b) 



In the following cases the return probability behaves 
as situations already investigated: 



(a) When £'^ = (i.e., 
inSec.IVC for Bn 



1] — 0) the results are presented 
= 0. 



(b) When Es » Ago (i.e., ry — >■ n/2) the behavior is 
as the situation presented in Sec. IV A for nonzero 
magnetic field. 

Away from these two limits, for any of the four (0,2) 
prepared states having one electron in each doublet (R, -I- 
and R, — ) the return probability goes smoothly from type 
8^ (at Es = 0) to type 8| (at E^ > A^o). The situation 
is not so interesting since the saturation value Poo is al- 
ways 0.50. The shape of P{ts) and the decaying times 
Td for each prepared state depend on the effective 8x8 
evolving (1,1) Hamiltonian. The result falls in between 
the two above mentioned limits (a) and (b). As shown in 
the previous sections the Hamiltonian for Es = is non- 
symmetric, while it is symmetric for Eg 3> Ago. In the 
latter case, P{ts) is independent of the prepared state 
and a smaller decaying time is observed. 
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FIG. 10. Case D. Results for a double dot with Ago 7^ as 
a function of the Zeeman energy produced by the perpendic- 
ular magnetic field, (a) Variances of the effective hyperfine 
field, (b) Saturation value of the return probability, Poo, and 
decaying time, Td, for the cases when the prepared state is 
the ground state or the highest energy excited state. 

On the other hand, the value of the saturation return 
probability changes if the (0,2) ground state or the (0,2) 
highest excited state is prepared. The ground state is the 
following Slater determinant: 



* 



(0,2) 



iv) 



R,K',- 
R,K,- 



1 

71 



rpLRcvalrjispin 



1 



+ ^cos?7|T^Rt™1s"p'") 
v2 



■ sin?7 



rpLRoval / rpspm 



rpS 



, (70) 



At zero-field rj — O and the prepared state is T'l^SKt.K'i, 
and we find type 4^ behavior with Poo = 1/3. For a 
dominant Zeeman energy 77 — 7r/2 and the prepared state 
is the spin polarized valley singlet, t^^S^'TIp'", with 
the spin triplet along the cc-direction; we then find type 
45^ behavior with Po^ =3/8. 

For the intermediate magnetic field regime, with the 
(0,2) state having two electrons in one doublet, the prob- 
lem is mapped to a double dot without the valley de- 
gree of freedom in a hyperfine field with the variances of 
Eq.(69b). Following Appendix A, we obtain the return 
probability by computing the averages of Eq.(AlO) with 
the probability distribution of Eq.(All) providing that 
CTh..^ ^Ch^.; =f?|| and o-^^ ,, =(T_L. 

We plot the standard deviations a± and crj| in Fig. 10(a) 
as a function of -Eg /Ago- There is an overall reduction of 
the HFI when the magnetic field increases, which explains 
the larger decaying time r^ plotted in Fig. 10(b). More- 
over, the standard deviation of the perpendicular com- 
ponent goes to zero, and in this limit the HFI becomes 
2-diniensional, approaching type 4" behavior. Thus, the 
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initial state in Eq.(70) becomes a valley singlet and the 
hyperfine interaction is unable to mix directly with the 
T™' partner. The saturation value of the return proba- 
bility is presented in Fig. 10(b). Poo can take any value 
between 1/3 and 3/8 as a function of the magnetic field. 
This interesting result allows a direct measurement of 
the spin-orbit coupling splitting and the hyperfine inter- 
valley mixing (see x and j/-components of the effective 
HFI in Eq.(68)) and its relation with the valley conserv- 
ing hyperfine mixing (z-component of the effective HFI). 



CONCLUSIONS 



TABLE III. Summary of the return probability results for the 
nine situations investigated. 



Poo 


Typo 


Cases 


n,. 


"(0,2) 


rdO-h 
for S-'T^p" 


/h 

for others 


1/3 


4= 


C 


4 


1 


- 


0.18935 


3/8 


4V 


A and C, B^O 


4 


1 


0.3218 


- 


0.40 


166 


A and B, B = 


16 


6 


0.185 


0.149 


1/2 


2= 


B and C, B^^iO 


2 


1 


- 


0.2944 


0.5 


8| 


^,^^0 


8 


4 


- 


0.238 


0.5 


8S 


C, B = 


8 


4 


0.3014 


0.2011 


0.5 


42 


C, B = B4 


4 


2 


- 


0.2585 


0.5 


4? 


C. B = B8 


4 


2 


0.4366 


0.2327 


1 


2? 


B and C, B^^iO 


2 


1 


oo 


- 



We have analyzed the expected return probabilities 
for a dephasing measurement in clean carbon nanotube 
based double quantum dots. We have focused on the 
intrinsic properties and therefore neglected disorder in- 
duced valley mixing and also Coulomb exchange, which 
are predicted^^ and measured^ to be small in multi- 
electron dots. In a forthcoming publication, we study 
the influence of valley mixing. 

We have shown that a multiple number of scenarios 
exists for the return probability experiment, due to the 
valley degree of freedom (as in Si-based DQDs,^^) which 
makes the system very different from a double dot in a 2- 
dimensional electron gas (2DEG). Here, more specifically, 
these scenarios are due to: (i) the non-trivial structure 
of the hyperfine coupling with the ^'^C nuclei that affects 
both the electron spin and valley degrees of freedom; (ii) 
the experimental preparation protocol that determines 
which of the six (0,2) states is prepared; (iii) the avail- 
ability of sixteen (1,1) states for the system in the evo- 
lution stage; (iv) the change (for every possible prepared 
state) of the subset of (1,1) states accessible in the evo- 
lution stage, and (v) the manifold of six possible return 
(0,2) states. The last point is an important difference 
to the 2DEG-bascd double dots, where only spin singlet 
returns to (0,2). Here the projection onto (0,2) is more 
generally determined by symmetry of the wavefunction, 
allowing only even left-right components to return. The 
level structure of the sixteen (1,1) and six (0,2) states de- 
pends on the values of the spin-orbit coupling and of the 
external magnetic field, through the Zecman interaction, 
diamagnetic effects, or both. 

In a 2DEG-based double dot the return probability 
shows two different behaviors being, type 4^ dynamics 
(Poo = 1/3) for zero field, or, type 2\ dynamics (Poo = 1/2) 
in the high magnetic field limit (we use the labeling intro- 
duced in Sec. IV). Here, depending on the parameters we 
find seven additional types of dynamics leading to satu- 
ration values Poo =3/8, 0.4, 1/2 and 1. The results for 
all the nine situations are presented in Table III. 

Type IGg dynamics can be found for zero magnetic 
field in the absence of spin-orbit coupling. In the cases 
with Ago = 0, the behavior of a Zeeman interaction only 
system (type 8| and 4") is very different from the sit- 
uation with both Zeeman and diamagnetic effects (type 



2f and 2"). For nonzero spin-orbit, the breaking of the 
spin degeneracy replaces, for zero magnetic field, type 
165 behavior with type 84 and the well known type 4^. 
At finite magnetic fields (when considering both Zeeman 
and diamagnetic effects) situations 2|, 2", 4|, 4^, 4" and 
once again A\ can be obtained depending on the value of 
the magnetic field and on the prepared state. 

In only two out of these seven novel situations (types 
2" and 4"), the return probability is associated with the 
system returning to the original prepared (0,2) state (see 
point (v) above). In all the remaining cases the sys- 
tem can be measured and also prepared in more than 
one (0,2) state, and therefore the functional dependence 
return of the probability on Ts depends on both the 
prepared state and the dynamics type. We have de- 
fined a shape-independent decaying time and we find, 
Td si0.149fi/ (T/i, in the fastest case (for zero spin prepared 
states in type IGg dynamics) and a infinite decaying time 
(since the system does not decay for type 2" dynamics) 
in the slowest case. It should be noted that we have as- 
sumed throughout that the time scale of the experiment 
is much smaller than Ti, the inelastic dephasing time; 
otherwise inelastic processes would relax the system to 
the ground state invalidating the investigation of the de- 
phasing in the return probability experiment as purely 
due to the hyperfine interaction. Therefore, the transient 
that defines the decaying time presented in Table III, t^, 
is to be understood as valid only for evolving times, r^, 
smaller than Ti. 

In addition to the robustness of case 2" we find, for 
types IGg, 84, 4" and A2, asymmetries and/or long de- 
caying times if the prepared state is a spin polarized val- 
ley singlet. The reason is that the hyperfine coupling 
of Eq.(36) does not introduces direct matrix elements 
between the LR-even spin polarized valley singlets and 
their LR-odd partners (see Eq.(53)); these two states can 
only be mixed by an inhomogeneous (in LR space) time- 
reversal symmetric term (i.e., spin-orbit coupling like), 
which does not appear in the HFI. In Rcf. 26 it has been 
shown that this intrinsically anisotropic hyperfine field 
gives rise to a dip in the spin-blockade signal as a func- 
tion of the orbital field. 

Here, we have shown that this property of the hyperfine 
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field also leads to an interesting behavior of the return 
probability when the (0,2) ground (or the highest excited) 
state is prepared for the case with nonzero spin-orbit cou- 
pling and the magnetic field is perpendicular to the tube 
axis, B±. As a function of the Zeeman energy, Eg, the 
groundstate changes from a spin-unpolarized state (for 
Es = 0) to a spin-polarized valley singlet (for Eg/Ago ^ 1 
) and the saturation return probability goes from 1/3 
(type 4^, i.e., effective 3-dimensional HFI) to 3/8 (type 
4", i.e., effective 2-dimensional HFI in the valley double 
dot). Measurement of Pqo and t^ as a function of B± 
would test the validity of the hypcrfinc Hamiltonian in 
Eq.(36) allowing, in principle, for the determination of 
the spin-orbit coupling and the hyperfine strength a^- 

Only a single return probability experiment^ has been 
reported in a carbon-based double dot. The result, only 
available for zero magnetic field, was an unexpected small 
return probability Poo ~ 1/6, that cannot be explained 
within the model presented here. We have shown that the 
minimum saturation return probability for coherent mix- 
ing is 1/3, similar to the situation in a spin-only double 
dot. Incoherent mixing will also not explain the experi- 
mental findings, since there the minimum return proba- 
bility is 1/4, which could happen for crossings type 4" or 
4\. We also note that by having worked in the high de- 
tuning limit in which the tunneling exchange is negligible 
we have obtained lower bounds of Poo , since it is known 
that this coupling reduces the effectiveness of the hyper- 
fine mixing and thus increases Poo-^'^ One could speculate 
that valley mixing is responsible for the discrepancy. In a 
forthcoming publication, we discuss the role of such mix- 
ing, which however also cannot explain the small ratio 
between P(0) and Poo seen in experiment. 

Clearly more experimental work is needed to better 
understand the rather rich structure of the carbon based 
double dots system, including the dependence of P{ts) 
on magnetic field. One interesting aspect would be to 
design alternative preparation protocols for being able 
to select different initial (0,2) states. 
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the simple Hamiltonian: 



He' = 



2Sh 



(0) 



2Sh 



(0) 







(Al) 



This Hamiltonian is valid for four out of the six Tq^, 
S'^^ pairs, specifically the zero-spin cases (see Eq.(49) 
and Eq.(53)). Here we use the notation of Eq.(25) and 
Eq.(20) for two single-particle single dot eigenstates with 
quantum numbers n and n' {n and n' must correspond 
to solutions with opposite spin projections). 



Following Eq.(28) for t^ = the state 
tialized. The system evolves as. 



LR„ 



V 



IS im- 



|e(r.)) 



'2(5/il°V, 



cos 



LR„ 



2Sh 



(0). 



lr: 



dd 



(A2) 



Therefore, for the realization r^ of the hyperfine field, 
the probability to find the system in the LR-even com- 
bination (i.e., to measure an (0,2) charge state after the 
adiabatical joining stage) is just: 



Prl{Ts) 



2Sh 



(0), 



(A3) 



We have to average the last oscillating function over the 
normal distribution that describes the hyperfine field. 
The inhomogeneous HFI component is given by 6hz — 
[h], i^ — h\ l^)/2 and the standard deviations for the com- 
ponents h], L and h], ^ are both \/2ah ■ Then the standard 
deviation of the Gaussian distribution for the frequency 



variable uj = 25h 



P'Hts) 



(0) 



/h is (T, 
1 



v27r(T,. 
1 + e-2(-< 



= 2(Jh/h. The final result is, 
duo cos^ [i^Tg) e ^"^'^ 



-fi(GhTs/hY 



■ (A4) 



We acknowledge useful discussions with S. Weiss, K. 
Grove-Rasmussen, H. Churchill, F. Kuemmeth, M. Lei- 
jnse, C. Marcus, B. Trauzettel and G. Burkard. 



2. Mixing of a Tq^ state with three S^^ states 
Analytical approach 



Appendix A: Analytical calculation of P{ts 



1. Mixing of a Tf,"", S^^ pair 



This case (type 2\ in Sec. IV) is valid whenever a LR- 
symmetric state and its partner LR-antisymmetric are 
mixed by the hyperfine interaction and no other states 
are involved the evolution Hamiltonian H^- In such a case 
the dynamics in the evolution subspace Sc is governed by 



The effective HFI Hamiltonians He' and He' pre- 
sented above (and also the intermediate situations we 
find for Ago in a perpendicular magnetic field) can be 
mapped to the problem of dephasing in a non-valley de- 
generated DQD as the one given in Eq.(37) and Eq.(39c). 
Here we present the derivation of the latter case and then 
we particularize for the three mentioned cases. 

The electron spin in each dot (^ = L, R) follows the 
evolution operators. 



U (ts) = cos (uJ^Ts) (To — i <T • n^ sin {uj^Ts 



(A5) 
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that describe precession around the direction of the hy- 
perfine field with frequencies, 



^« = T\Jihx,if + {hy^(f + {h^.cf 



(A6) 



The normahzed vectors in Eq.(A5) point in the direction 
of the local hyperfine field. 



n? 



hwf 



{hx.Cihy^^, h 



z,a 



(A7) 



Using the former evolution operators it follows that an 
(1,1) Slater determinant in the double dot evolves as. 



U{ts) 






cricr2 



K,aArs)U^.,,{Ts) 



R,(T2 

L,cri 



(A8) 



At Ts—0 the system is initialized in Tq^S^p'", the only 
available LR-symmetrical state in Eq.(38). In order to 
time evolve the last two-particle state, we use its Slater 



determinant version, ( -A ) 



^\\) /V2. We apply the 



evolution operator of Eq.(A8) and project the result back 
to the initial LR-even state. 



V^J^) = (T^i^S^P^-l U{r,) iT^i^S^P'") = I {U^^Ul, + Uf,U^^ - U^.Uj;^ + Uf^U. 



(A9) 



cos(a;RTs)cos(wLTs) 
sin(u;LTs) sin(a;Rrs) 



{hx.nhx.L + hy,j^hy^i, + hz^j^hz.h) ■ 



The probability of finding the system in the original state — i.e., of measuring its (0,2) partner state after the adiabatical 
joining stage — is then the square of the former amplitude. 

Since the probability density functions of the hyperfine field components (Gaussian distributions with zero mean) 
are even, the odd powers terms in those components within PniTs) do not contribute to the average. We arrive to 
the well known expression: 

(|f)'sin2(c.Lr.)^ 

3=x,y,z L \ ^ '"^ / L-HF \ ^ "' / R-HFJ 

(AlO) 
where (Q (/ix,{, ^y,{j ^2,^))c_jjp stands for the average of the function Q {hx,^, hy^^, hz.^) over the hyperfine fields of 
the ^ (L or R) dot. The probability density function (in the 3-dimensional space of hx^^,hy,^ and hx.^) is 



P{ts) ^ (cos^ (a;LT^))L_HF (*=°^^ ('^R^''))r-hf+ E 



— ) sin^ (ujrTs) 
wr/ , 



FsD {hx,^, hy^^,hz^s) 



1 



1 



1 



exp 



1 / 


f Ki ' 


2 

) + 


^ hy,i ' 


2 

+ 


^ Ki ' 


)11 


A 


K'^K.^^, 


/ 


V^'»„,«/ 




y^h^.i^ 


' l\ 



(All) 



Ji-KOh^i, v27r(Tft^_5 V27rcr^^^ 

Where we have used the following standard deviations of the Gaussian distributions for the HFI components, 

<^h.^,i = Std(/i2:,5) , (Thy^^ ^ SiA{hy^^) , o-ft^_j = Std{hz,{) ■ (A12) 

We see below that the degree of anisotropy arising from a difference in the last quantities affects the averages in 
Eq.(AlO) and therefore the return probability. 



a. Statistical isotropic 3-dimensional effective hyperfine field 



4 , , F 

When the effective Hamiltonian is Hc^ the three hyperfine components share the same standard deviation, a\s 



CTh^ 5 =(7hy^^ =<^h,^^ and therefore, 

/OO nOO nOO 

-co J —QO J —CO 



OO /'OO 

OO 
TT />27T /'+00 

d0 d^ drQ (r, 9, ^) Fgo (r, 6, ip) , 



^0 

2 , 







, ^ , r^ sin(6') 
FsD (r, e, if) = -^ exp 



(72^4)' 



2 (T 



(A13) 
(A14) 



where we have made a change to spherical coordinates with r'^ = {hx,^) + (hy.^) + {h^^^Y 
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The two types of averages that appear in Eq.(AlO) are 
obtained by integrating, 



C3d(ts) = (cos^ {uj^Ts)) 



C-HF 



l(i + 9(Ts4jh 



S3d{Ts) = 




sin^ (uj^Ts) 



1 - g (tsO-I^/H 



-HF 



]{x) = e-2^' (l - 4x2) 



(A15) 



(A16) 
(A17) 



Since hu^ = r the arguments of the sinusoidal func- 
tions in spherical coordinates is rTs/h. Note that in 
SsD the average is independent of the direction of the 
hyperfine component {j = x,y,z). This is valid here be- 
cause the effective hyperfine field is statistically isotropic. 
Then, for simplicity, the integral is computed using the 
^-component, hj^^ = r cos{9) . 



In the investigated situations in Sec. IV the standard 
deviations are equal in the two dots, a^ =a^B =(t^s ; the 
return probability becomes 

P4i (r,) = ^ (1 + 5 (rs^ii/h) + g^ (rsa^/h)) . (A18) 



b. Statistical isotropic 2- dimensional effective hyperfine 
field 

As discussed in Sec. II A 1 the effective Hamiltonian 
Hc^ can be mapped to the GaAs zero- field double dot 
but it must be assumed that the effective 2;-component 
of the hyperfine field is absent. As h^^^ is identically 
zero we must not average over it, therefore, instead of 
the probability density function given in Eq.(All) a two 
dimensional probability density function must be used. 
Adding the fact that the standard deviations of the in- 
plane components are identical, crfn — cr;;^^ = '^hy^^^ it 
becomes useful to work in polar coordinates. The aver- 
ages are obtained as follows. 



(Q(/lx.Ci^2/,«))^-HF 



27r /' + 00 

df / drQ (r, ip) F2D {r, p) , 



F2D {r, (/?) 











V2t 



exp 



2 C7 



(A19) 
(A20) 



Then we define the 2-dimcnsional averages needed to 
compute the return probability as: 



'hi 



S2d{Ts) 



hw 



sin (uj^Ts) 



(A21) 
3 = X, y. 



I 

As in 5*30 the average in ^20 is independent of the di- 
rection of the hyperfine component (j = x,y). This is 
valid here because the effective in-plane hyperfine field 
is statistically isotropic. Then, for simplicity, the inte- 
gral is computed using the x-component, hj^^=rcos{if). 
The obtained results are presented and discussed in 
Sec.IIAl. 
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